arXiv:gr-qc/9706083v3 18 Dec 1998 


PHYSICALLY VALID BLACK HOLE INTERIOR MODELS 


Giulio Magli 


Dipartimento di Matematica del Politecnico di Milano 
P.le Leonardo da Vinci 32, 20133 Milano, Italy 
tel: 39-2-23994538 (voice) 39-2-23994568 (fax) 
e-mail magli@mate.polimi.it 


New, simple models of “black hole interiors”, namely spherically symmetric solutions 
of the Einstein held equations in matter matching the Schwarzschild vacuum at spacelike 
hypersurfaces “R < 2M” are constructed. The models satisfy the weak energy condition 
and their matter content is specihed by an equation of state of the elastic type. 
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Spherically symmetric, static models of stars in General Relativity are composed by a 
regnlar solntion of the Einstein held eqnations in matter matching smoothly the Schwarz- 
schild vacnnm at a snrface which lies ontside the event horizon, while the Schwarzschild 
black hole is nsnally regarded per se as a singnlar solntion of the Einstein held eqnations; 
its physical acceptability from the astrophysical point of view relies in the fact that the 
singnlarity is covered by the horizon. However, of conrse, the problem of singnlarities is 
far from being solved, and it is generally believed that qnantnm ehects shonld restrict 
the cnrvatnre to hnite valnes. As a conseqnence, in recent years the possibility of an 
alternative pictnre in which also black holes have non-singnlar interiors has been actively 
investigated [1-6]. For “interior strnctnre” of a (non-rotating) black hole one means a 
solntion of the Einstein held eqnations in matter matching smoothly the Schwarzschild 
solntion in a region which lies inside the Schwarzschild radins. Of conrse, specnlations of 
this kind are a priori academic becanse no one can commnnicate to an external observer 
the discovery of a “closed world” inside a black hole. We consider, however, interesting to 
check whether General Relativity allows the existence of this kind of objects at a pnrely 
classical level. 

Investigating the possible “interior strnctnre” of black holes, one is immediately faced 
with rigorons theorems [7-9] which impose severe restrictions: no solntions satisfying the 
strong or the dominant energy condition may exist. Therefore, the physical acceptability 
of snch models mnst be restricted to the weak energy condition only. Models of regnlar 
black holes interiors satisfying this condition have been fonnd by Mars, Martin-Prats and 
Senovilla [6]. Recently, Horde [10] remarked that a model of this kind is already known 
from nnpnblished work by J. Bardeen (see also [11]). Borde has shown that the mechanism 
responsible for avoidance of singnlarities in all regnlar black holes is topology change. 

Taken together, the above mentioned resnlts give ns a qnite complete nnderstanding 
of “regnlar black hole physics”. It is, however, worth mentioning that, if we really want 
to check if natnre allows the existence of snch objects, we shonld hnd models satisfying 
an eqnation of state, so that a clear physical interpretation of the matter content of their 
sonrce can be given. As far as the present anthor is aware, no examples of this kind are 
known, and onr aim here is to hll this gap. Onr starting point is that strongly collapsed, 
relativistic bodies snch as nentron stars typically exhibit elastic-solid properties dne to a 
process of crystallization of matter at very high densities (see e.g. [12]). Therefore, there 
is a possible interpretation of isolated, anisotropic sonrces of strong gravitational helds as 
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highly collapsed, elastic solid objects. In the present paper, we propose a very simple class 
of “elastic-solid” black hole interiors. 

Following Ref. [6], we consider the spherically symmetric line element in Eddington- 
Finkelstein coordinates: 


ds^ = —e'^^xdu^ + 2e^^(i'U(ir + r^{d9^ + sm‘^{9)d(l)^) , (1) 


where (3 and x are fnnctions of r only. It is nsefnl to represent x as 


X{r) 


2^(r) 


The form (1) of the line element is very convenient in dealing with the black hole interior 
problem. In fact, observe that one can always nse the transformation between the nnll 
coordinate u and the Schwarzschild “time” t given by 

dt{u, r) = du - -—dr 

X{r) 

which brings the line element in its Schwarzschild form, bnt this transformation is allowed 
only separately within regions in which the fnnction x has a constant sign. It is worthwhile 
to note that, rewriting (1) as 


ds‘^ — —e^^X 




+ —dr^ + r‘^{d9‘^ + sin^{9)d(l)^) , 
X 


( 2 ) 


the line element in Schwarzschild form is static if x is positive, while it is non-static if x 
is negative. We shall refer towards to snch regions as to R— and T— regions respectively 
[13]. The line element (1) can be nsed to describe spacetimes which contain both R— and 
T— regions. 

To proceed fnrther, we assnme as a sonrce an isotropic elastic medinm. The energy- 
momentnm tensor of an elastic-solid sphere has been thoronghly discnssed in Ref. [14]. In 
particnlar, it was shown that a spherically symmetric, relativistic body can be characterized 
by three qnantities e,p, O. In terms of snch qnantities the eigenvalnes of the energy- 
momentnm tensor may be written as 


Ao — —e , 

Ai = p + 2hl , 

A2 = As = p — O . 


( 3 ) 
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In i?—regions, Aq corresponds to the time-like eigenvector and therefore e is the energy 
density, while Ai is the radial stress, which the snm of two terms: the isotropic part p and 
the qnantity Q which is a measnre of the response of the body to strains that change its 
shape withont changing is volnme: if O vanishes, the elastic medinm is simply a perfect 
flnid. The eigenvalne A 2 always corresponds to a space-like eigenvector, and therefore is 
the tangential stress, while, in T—regions, the role of Aq and Ai is reversed. 

The Einstein held eqnations may be written as follows: 

p = dyrr^e , 

X/3'= 27rr(e + p + 20) , 

2 x(^p' + 20 ' + 6^^ =-(e + p + 20)(4x/3' + x') , 

where a dash denotes the derivative with respect to r. This system is composed by three 
eqnations for the hve nnknowns (x, P, e,P, hi) and therefore it is not closed nntil we snpple- 
ment it with an eqnation of state. Generally speaking, to specify the eqnation of state for 
an isotropic elastic body one gives the energy density as a fnnction of the three invariants 
of the strain tensor. The stresses may then be calcnlated as derivatives of the energy with 
respect to snch invariants (in the case of a perfect hnid, this obvionsly amonnts to say 
that the eqnation of state relates energy and specihc volnme n, and that the pressnre may 
be obtained as the derivative of the energy with respect to v). In the particnlar case of 
spherical symmetry, it has been shown [14] that assigning the state eqnation of the body 
consists in giving the energy density as a fnnction of only two invariants of the strain 
tensor, one of them being the specihc volnme (so that the derivative of e with respect to 
it governs the isotropic response p) and the other one being the qnadratic invariant of the 
deformation, the “response” associated to this invariant being O. 

The weak energy condition reqnires 

> 0 , (5) 

for any timelike vector v^. Taking into acco un t (3), the above ineqnality leads to 

e>0, e + p + 2O>0, e + p — O>0, 


p T 20 ^ 0 , e T p T 20 ^ 0,0^0, 


in i?—regions, and to 
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in T—regions. 

A very simple way to simplify the above conditions is to consider that particnlar class 
of elastic materials satisfying 

e + p + 2n = 0 . (6) 

In fact, the weak energy condition for snch materials rednces to 


e > 0 , O < 0 , 


(7) 


in any region. We stress that the above choice, besides being very convenient from the 
mathematical point of view, has a clear physical meaning. It is, in fact, eqnivalent to a 
“constitntive” partial differential eqnation for e. The solntion of snch eqnation identihes a 
well-dehned class of elastic materials which depends on the choice of an arbitrary fnnction 
of one variable only (we refer the reader to Ref. [14] for details). This arbitrariness may 
be nsed to assign the “on shell” valne e = e(r). The system (4) then becomes closed and 
rednces to 

p = dyrr^e , 


/ 9 '- 0 , 



Therefore, the fnnction y is given by 


( 8 ) 


y(r) 1 - — f s‘^e{s)ds , 

r Jo 


the fnnction /? is a constant which may be set eqnal to zero withont loss of generality, and 
the line element assnmes the simple form 


= —x{r)du^ + 2dudr + r‘^{d9‘^ + siv?{9)d(f?) . 


(9) 


This solntion may be seen as a “variation of the mass” of the Schwarzschild solntion and 
has many interesting featnres. For example, it is a Kerr-Schild geometry and it arises 
natnrally as the spherically symmetric limit of a recently fonnd class of sonrces of the 
Kerr metric [15]. It arises also in completely different contexts, for example in general- 
relativistic electron models (see e.g. [16] and references therein) or in vacnnm polarization 
[5]. 

Let ns now consider the conditions which have to be imposed on the free f un ction e 
in order to obtain a regnlar black hole model. 
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First of all, the weak energy condition (7) reqnires e to be a non-negative fnnction, 
while O = re'/6 mnst be negative, so that e mnst be a decreasing fnnction of r. The 
regnlarity of the spacetime at r = 0 reqnires ~ as r tends to zero, and therefore e mnst 
be hnite there. Finally, in order to obtain a black-hole model, we want to match smoothly 
the interior solntion (9) with the Schwarzschild solntion at a space-like hypersnrface r = R. 
The matching between two metrics is smooth (no snrface distribntions of matter-energy 
arise) if the hrst and the second fnndamental form are continnons at the matching snrface. 
Continnity of the metric reqnires the mass of the vacnnm solntion to be given by 

M := pi{R) . 


To obtain a black-hole model, we reqnire R to belong to a T—region: 


2M 


< 0 . 


( 10 ) 


Continnity of the second fnndamental form reqnires vanishing of the energy density since 
r = R is a space-like hypersnrface. The energy density is — (p + 20) in a T—region, bnt 
hnally in onr model this eqnals e, so that we reqnire e{R) = 0. 

Regarding the strnctnre of the solntions, it is easy to check that the fnnction y mnst 
vanish at at least one valne Tq of r, becanse x(r) is continnons and satishes x(0) = 1, 
x(R) < 0. This is a characteristic shared by all regnlar black holes interiors satisfying the 
weak energy condition [6]. A nice featnre of onr solntions is that there exist always only 
one valne Tq, as a graphic comparison between r and the monotone increasing fnnction 
2|u(r) immediately shows. A simple calcnlation also shows that the Riemann invariants are 
hnite, and therefore the solntion is non-singnlar and Tq is a Canchy horizon. At the horizon, 
the Plebahski type of the energy-momentnm tensor degenerates from [25'i — S 2 — T] to 
[25” — 2A^] (the timelike and the radial spacelike eigenvector degenerate into a donble nnll 
eigenvector). It may, however, be easily shown that the weak energy condition remains 
satished there. 

The maximal analytic extension of onr solntions is similar to that already well known 
for others “black hole interiors” (see e.g. [5-6]). It closely resembles that of the Reissner- 
Nordstrom spacetime with the key difference that the singnlarity at r = 0 is replaced by 
the matter-hlled region. 

It may be instrnctive to discnss some explicit examples. Considering a power-law 
distribntion of density 



e = a 


( 11 ) 
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(where a and n are positive quantities), the Schwarzschild mass is 

4n 


M = 


-ttR a , 


3(n + 3) 

and the solution is a black hole interior provided that 

■3(n + 3)^^/^ 


R > -Rmin ■ 


Snna 


therefore, at hxed n, the minimal required value of r at the matching hypersurface decreases 
as a increases. The solution is given by (9) with 


X{r) 


1 - 


2M / r \ 2 r 


nR \R 


n + 3 — 3 


R 


Since the central value a = e(0) is really the density at the centre because r = 0 
always belongs to a i?—region, it seems reasonable to assume it around the typical values 
for neutron stars, namely of order 10^® g/cm^. Considering, for instance, the case of a 
linear function e (n = 1 in (11)) we can calculate Since we are using relativistic units, 

we can express it in Kilometers (it is about 25 Km), although it would be more correct 
to express it in time-units (about 8.2 x 10“^ seconds) because it is in fact an instant of 
time, belonging to a T—region. In any case, of course, the unique quantity which is really 
measurable by an external observer is the mass of the Schwarzschild black-hole. Assuming 
R to be 26 Km one obtains this mass to be M ^ 8.5Mq, with a Schwarzschild radius 
2M Ri 28 Km. 


The author thanks Jerzy Kijowski for useful discussions and comments. 
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